* formal_graph.do	JEP	01/27/15
* Build regression table
* Formal sector, wages

*************************** Change log ********************************** 

/*
	
*/	

clear 
version 13.1

* Project info
cap project, doinfo
if _rc==198 {
	if c(username)=="J16339" loc master "B:\Col_Minimum_Wage"
	else if c(username)=="A15858" loc master "B:\Col_Minimum_Wage"	
	else if c(os)=="Windows" loc master "C:/Users/jorpp/Dropbox (Brown)/Col_Minimum_Wage"
	else loc master "/home/jperez/AA_Minimum_Wage"
	loc pr=0
}
else {
	local master "`r(pdir)'"
	local doname "`r(dofile)'" 
	* Project calls
	project, relies_on("`master'/Code/Est/rifregtable2.do")
	project, relies_on("`master'/Code/Est/graph_rif.do")
	foreach x in 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 {
		project, uses("`master'/Ster/rif4_new_indcc/reg_emp_bp_tq(1999q1)_bite_salario_mensual_real_50_edom_indcc_p`x'.ster")
	}
	loc pr=1
}

cap program drop mex
program define mex
	cap mkdir "`2'/Doc"
	cap mkdir "`2'/Doc/tosubmit2"
	cap mkdir "`2'/Doc/tosubmit2/fig"
	graph export "`2'/Doc/tosubmit2/fig/`1'.pdf", replace
	graph export "`2'/Doc/tosubmit2/fig/`1'.png", replace
	graph export "`2'/Doc/tosubmit2/ILR-Figs/Fig4_`1'.png", replace	
end	

do "`master'/Code/Est/graph_rif.do"

*graph drop _all
glo qtiles "5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90"
set scheme s1color
loc i=1
foreach spec in /*reg_emp_bp reg_emp_bq */ reg_emp_bp {
	if `i'!=1 loc leg "legend(off)"
	if "`spec'"=="reg_emp_bp" loc subt "Bartik price shocks"
	if "`spec'"=="reg_emp_bq" loc subt "Bartik quantity shocks"
	if "`spec'"=="reg_sst_emp_bp" loc subt "Bartik price shocks - City trends"
	if "`spec'"=="reg_emp_bp_ixc" loc subt "Bartik price shocks - City x Industry FE"
	if "`spec'"=="reg_sst_emp_bq" loc subt "Bartik quantity shocks - City trends"
	if "`spec'"=="reg_emp_bq_ixc" loc subt "Bartik quantity shocks - City x Industry FE"
	
	else loc leg "legend(label(2 "UQE") label(1 "95 % CI"))"
	graphrif, master("`master'") ster(Ster) folder(rif4_new_indcc) qtiles($qtiles) spec(`spec') time(tq(1999q1)) rareaopts(fcolor(gs13) lcolor(none) lpattern(dash) ) twowayopts( xtitle(Percentile) ytitle(Coefficient) yscale(range(-0.01 0.01)) xscale(range(5 90)) ylabel(-0.01(0.005)0.01, angle(horizontal)) legend(off) name(gf`i'), (pci -0.0125 23 0.01 23, lcolor(black)) (pci -0.0125 31 0.01 31, lcolor(cranberry) lpattern(dash)))   lineopts(lcolor(black) xlab(10(10)90)) 
	loc i=`i'+1
}

* Had to add text by hand
graph play "`master'/Code/Doc/tosubmit2/formal_edit.grec"

mex formal "`master'"

if `pr' project, creates("`master'/Doc/tosubmit2/fig/formal.pdf")

/*
do "`master'/Code/Est/graph_rif_alternativo.do"
foreach spec in /*reg_emp_bp reg_emp_bq */ reg_emp_bp {
	if `i'!=1 loc leg "legend(off)"
	if "`spec'"=="reg_emp_bp" loc subt "Bartik price shocks"
	if "`spec'"=="reg_emp_bq" loc subt "Bartik quantity shocks"
	if "`spec'"=="reg_sst_emp_bp" loc subt "Bartik price shocks - City trends"
	if "`spec'"=="reg_emp_bp_ixc" loc subt "Bartik price shocks - City x Industry FE"
	if "`spec'"=="reg_sst_emp_bq" loc subt "Bartik quantity shocks - City trends"
	if "`spec'"=="reg_emp_bq_ixc" loc subt "Bartik quantity shocks - City x Industry FE"
	
	else loc leg "legend(label(2 "UQE") label(1 "95 % CI"))"
	graphrif, master("`master'") ster(Ster) folder(rif4_new_indcc) qtiles($qtiles) spec(`spec') time(tq(1999q1)) rareaopts(lcolor(black) lpattern(dash) ) twowayopts( title("Formal") xtitle("Percentile of the wage distribution") ytitle("Effect on wages for each p.p. of vinculation") yscale(range(-0.01 0.01)) xscale(range(5 90)) ylabel(-0.01(0.005)0.01, angle(horizontal)) legend(off) name(gf2`i'), (pci -0.0125 23 0.01 23, lcolor(black)) (pci -0.0125 31 0.01 31, lcolor(cranberry) lpattern(dash)))    lineopts(lcolor(black) xlab(10(10)90)) 
	loc i=`i'+1
}

graph play "`master'/Code/Doc/tosubmit2/formal_edit.grec"

*/

graph rename formal
